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Oh Abstract. In most cosmological models, primordial black holes (pbh) should have 

formed in the early Universe. Their Hawking evaporation into particles could even- 
I tually lead to the formation of antideuterium nuclei. This paper is devoted to a 

' first computation of this antideuteron flux. The production of these antinuclei is 



studied with a simple coalescence scheme, and their propagation in the Galaxy is 
treated with a well-constrained diffusion model. We compare the resulting primary 
flux to the secondary background, due to the spallation of protons on the interstel- 
lar matter. Antideuterons are shown to be a very sensitive probe for primordial 
black holes in our Galaxy. The next generation of experiments should allow in- 
vestigators to significantly improve the current upper limit, nor even provide the 
first evidence of the existence of evaporating black holes. 

Key words. Black hole Physics - Cosmology: miscellaneous 

1. Introduction 

Very small black holes could have formed in the early Universe from initial density 
inhomogeneities fHawking I1971f) . from phase transition (Hawking I1982fl . from collapse 
of cosmic strings (Hawking I1989|l or as a result of a softening of the equation of state 
fCanuto I1978|l . It was also shown by Choptuik fChoptuik I1993[) and, more recently, 
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studied in the framework of double inflation fKim |2"UUU|) . that pbhs could even have 
formed by near-critical collapse in the expanding Universe. 

The interest in primordial black holes has been revived in the last years for several 
reasons. On the one hand, new experimental data on gamma-rays ( Connaughton Il998|l 
and cosmic rays (Barrauei al. I2()02[l together with the construction of neutrino detectors 
(Bugaev & Konishchev l2U(Jlf) . of extremely high-energy particle observatories (Barrau 
200dJ and of gravitational waves interferometers (Nakamura et al. I1997f) give interesting 
investigational means to look for indirect signatures of PBHs. On the other hand, pri- 
mordial black holes have been used to derive interesting limits on the scalar fluctuation 
spectrum on very small scales, extremely far from the range accessible to CMB studies 
(Kim et al. 119991 Blais et al. I2003JI . It was also found that PBHs are a useful probe of 
the early Universe with a varying gravitational constant (Carr|3p00). Finally, significant 
progress has been made in the understanding of the evaporation mechanism itself, both 
at usual energies (Parikh & Wilczek l2000(l and in the near-Planckian tail of the spectrum 
(Barrau & Alexevev IMH1 Alexeyev et aZ.ESUl Alexeyev et al. fM?^ . 

For the time being there is no evidence in experimental data in favour of the existence 
of pbhs in our Universe. Only upper limits on their number density or on their explosion 
rate have been obtained (Barrau et al. 120021 MacGibbon & Carr ll99l|l . As the spectra of 
gamma-rays, antiprotons and positrons can be well explained without any new physics 
input {e.g. pbhs or annihilating supcrsymmctric particles) there is no real hope for any 
detection in the forthcoming years using those cosmic-rays. The situation is very different 
with antideuterons which could be a powerful probe used to search for exotic objects, as 
the background is extremely low below a few GeV (Chardonnet et oZ. 119971 Donato et al. 
20003- Such light antinuclei could be the only way to find pbhs or to improve the current 
limits. This paper is organized along the same guidelines as our previous study on pbh 
antiprotons (Barrau et al. I2002f) . to which the reader is referred for a full description of 
the source and propagation model used. The main difference is the necessity to consider 
a coalescence scheme for the antidcuteron production. We compute the expected flux 
of antideuterons for a given distribution of pbhs in our Galaxy, propagate the resulting 
spectra in a refined astrophysical model whose parameters arc strongly constrained and, 
finally, give the possible experimental detection opportunities with the next generation 
of experiments as a function of the uncertainties on the model. 

2. Antideuterons emission 

2.1. Hawking process and subsequent fragmentation 

The Hawking black hole evaporation process can be intuitively understood as a quantum 
creation of particles from the vacuum by an external field. The basic characteristics can 
be easily seen through a simplified model (see Frolov & Novikov i|1998[l for more details) 
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which allowed Schwinger to derive, in 1951, the rate of particle production by a uniform 
electric field and remains correct, at the intuitive level, for black hole evaporation. If 
we focus on a static gravitational field, it should be taken into account that the energy 
of a particle can be written as E = —p^S,^, where p^ is the four- momentum and £ M is 
the Killing vector. The momentum being a future-directed timelike vector, the energy 
E is always positive in the regions where the Killing vector is also future-directed and 
timelike. If both particles were created in such a region, their total energy would not 
vanish and the process would, therefore, be forbidden by conservation of energy. As a 
result, a static gravitational field can create particles only in a region where the Killing 
vector is spacelike. Such a region lies inside the Killing horizon, i.e. the £ 2 = surface, 
which is the event horizon in a static spacetime. This basic argument shows that particle 
creation by a gravitational field in a static spacetime (this is also true in a stationary 
case) is possible only if it contains a black hole. Although very similar to the effect of 
particle creation by an electric field, the Hawking process has a fundamental difference: 
since the states of negative energy are confined inside the hole, only one of the created 
particles can appear outside and reach infinity. 

The accurate emission process, which mimics a Planck law, was derived by Hawking, 
using the usual quantum mechanical wave equation for a collapsing object with a post- 
collapsc classical curved metric instead of a precollapse Minkowski one (Hawking I1975f) . 
He found that the emission spectrum for particles of energy Q per unit of time t is, for 
each degree of freedom: 

d 2 ^ r s 



dQdt 



(e*p(Wl^)-(-l) 2s ) 



(1) 



where contributions of angular velocity and electric potential have been neglected 
since the black hole discharges and finishes its rotation much faster than it evaporates 
(MacGibbon & Webber 119751 Page 11977(1 . k is the surface gravity, s is the spin of the 
emitted species and T s is the absorption probability. If we introduce the Hawking tem- 
perature (one of the rare physical formula; using all the fundamental constants) defined 
by 

he 3 10 13 e 

T = — ^ « ^-^GeV (2) 

16irkGM M V ; 

the argument of the exponent becomes simply a function oiQ/kT. Although the absorp- 
tion probability is often approximated by its relativistic limit 
1087r 2 G 2 A/ 2 Q 2 

lim T s = -g-g 3 

Q^oo h^C 3 

we took into account in this work its real expression for non-relativistic particles: 
4^.(0 M, M ) a _ a 

where a s is the absorption cross-section computed numerically (Page 11976(1 and /j, the 
rest mass of the emitted particle. Even if this mass effect is partially compensated by the 
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pseudo-oscillating behaviour of the cross-section and remains at the level of a correction, 
we found some substantial discrepancies between the geometric limit and the numerical 
calculation which justifies this technical complication. 

As was shown by MacGibbon and Webber (MacGibbon & Webber Il990fl . when the 
black hole temperature is greater than the quantum chromodynamics confinement scale 
A-QCD, quark and gluon jets are emitted instead of composite hadrons. To evaluate the 
number of emitted antinuclcons ./V , one therefore needs to perform the following convo- 
lution: 

d 2 A^v V- f°° T sAQ,T) ( o_ ^.V 1 dg jR {Q,E) 



where otj is the number of degrees of freedom, E is the antinucleon energy and 
dgjjy(Q,E)/dE is the normalized differential fragmentation function, i.e. the number of 
antinuclcons between E and E + dE created by a parton jet of type j and energy Q. The 
fragmentation functions have been evaluated with the high-energy physics frequently- 
used event generator PYTHIA/JETSET (Tiostrand ll994[) . This program is based on the 
so-called string fragmentation model (developed by the Lund group) which is an explicit 
and detailed framework where the long-range confinement forces are allowed to distribute 
the energies and flavours of a parton configuration among a collection of primary hadrons. 
It has received many improvements related, e.g., to parton showers, hard processes, Higgs 
mechanisms and it is now in excellent agreement with experimental data. 

2.2. Coalescence scheme 

In the context of proton- nucleus collisions it was suggested that, independently of the 
details of the deuteron formation mechanism, the momentum distribution of deuterons 
should be proportional to the product of the proton and neutron momentum distributions 
(see Csernai & Kapusta (|1986|l for a review). This was based on phase space consider- 
ations alone: the deuteron density in momentum space is proportional to the product 
of the proton density and the probability of finding a neutron within a small sphere of 
radius po around the proton momentum. Thus: 
d 3 N d _ 4tt 3 / d 3 N p \ ( d 3 N,, 



where po is the coalescence momentum which must be determined from experiments. 
The very same arguments can be used for antideuterons resulting from an antiproton 
and antincutron momentum distribution. In our case, the coalescence scheme has to be 
implemented directly within the PBH jets as no nuclear collision is involved. We defined 
the following procedure: 

- for each hadronic jet resulting from a parton emitted by a pbh, we search for 
antiprotons 
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- if an antiproton is found within the jet, we search for antineutrons 

- if an antineutron is also found within the same jet, we compare their momenta 

- if the difference of the antiproton and antineutron momenta is smaller than the 
coalescence momentum po, we consider that an antideutcron should be formed. 

As the coalescence momentum po is not Lorentz invariant, the condition must be 
implemented in the correct frame, namely in the antiproton-antineutron center of mass 
frame instead of the laboratory one. Fig. Ogives the differential spectrum of antiprotons 
resulting from 1.9 x 10 8 u quark jets generated at 100, 75, 50, 25 GeV and the subsequent 
distribution of antideuterons with po = 160 MeV. The ratio between the antideuteron 
and antiproton spectra is of the order of a few times 10 -5 , which reflects the mean 
amplitude of the cosmic antideuteron flux from PBHs, given in Section 4 of this article, 
with respect to the one given for antiprotons in Barrau et al. (2002). This value is not 
surprising at it is in reasonable agreement with: 

- the Serpukhov experimental ratio of the p to D production cross-sections for proton- 
proton interactions measured at y/s = 11.5 GeV (between 1.9 x 10~ 5 and 3.5 x 10~ 5 
depending on the transverse momentum) as given in Abranov et al. (11987[) and for 
proton-aluminium interactions measured at the same energy (7 x 10~ 5 for a center of 
mass transverse momentum around 270 MeV) as given in Binon et al. (jEMO) 

- the CERN-ISR experimental ratio of the p to D production cross-sections for proton- 
proton interactions measured at \fs = 53 GeV (between 10~ 4 and 3.4 x 10 -4 ) as given 
in Alpcr et al. l|1973j> and Gibson et al. Ijl978jl 

- the theoretical cosmic secondary p to D ratio (around 10 -5 ) as given in Chardonnct 
et al. lfTM7|l 

- the theoretical cosmic primary p to D ratio from neutralinos (around a few times 
10" 4 ) as given in Donato et al. ( l2TjOTj|> 

- the simulated ratio of p to D fluxes created within the Earth's atmosphere (around 
10~ 5 ) as evaluated with the program (Derome et al. P000|) ) that was used to explain 
AMS experimental data (Derome, private communication). 

Although the orders of magnitude are correct, large discrepancies between these 
theoretical and experimental results can be noticed. This is taken into account in this 
work by allowing the coalescence momentum to vary between 60 MeV and 285 MeV, 
numbers than can be considered as "extreme" possible values. 



The flux of emitted antideuterons should now be written as: 



where dgjjy(Q 7 E,po)/dE is the fragmentation function into antideuterons evaluated with 
this coalescence model for a given momentum pq. As the mean number of produced 
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Fig. 1. Upper plot: antiproton differential spectrum obtained with 1.9 x 10 8 u quark 
jets generated by PYTHIA at 100, 75, 50, 25 GeV. Lower plot: antideuteron spectrum 
obtained in the same conditions with a coalescence momentum po = 160 MeV. 

antideuterons per jet is extremely low, millions of events were generated for each energy 
and each partonic degree of freedom. Some interpolations are also required to avoid a 
diverging computing time: the associated uncertainties have been found to be negligible. 



2.3. Convolution with the mass spectrum 

The above expression gives the antideuteron flux due to a single black hole of temperature 
T. As PBHs of different temperatures (or masses) should be present, this flux must be 
integrated over the full mass spectrum of pbhs: 
d 2 N B (M,E) d 2 n(r,z) 



q prtm (r,z,E) 



dEdt 



dMdV 



-dM 



with 

dn 

dM ' 



M 2 for M < AL 
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— ^ cx M~ 5/2 for M > M* 
dM 

where M* w 5 x 10 14 g is the initial mass of a PBH expiring nowadays. As explained in 
(Barrau et al. |2T)T)2| , the shape below M* does not depend on any assumption about 
the initial mass spectrum whereas the shape above Af* relies on the assumption of a 
scale-invariant power spectrum. The resulting distribution is, then, normalized to the 
local pbh density pq. The spatial dependence of this source term is given in eq. 

3. Propagation and source distribution 

The propagation of the antideuterons produced by pbhs in the Galaxy has been studied 
in the two zone diffusion model described in Donato et al. H'2()()1J) . 

In this model, the geometry of the Milky Way is a cylindrical box whose radial extension 
is R = 20 kpc from the galactic center, with a disk whose thickness is 2h = 200 pc and 
a diffusion halo whose extent is still subject to large uncertainties. 

The five parameters used in this model are: Kq, 5 (describing the diffusion coefficient 
K{E) = Kof3R d ), the halo half height L, the convective velocity V c and the Alfven velocity 
V a . They are varied within a given range determined by an exhaustive and systematic 
study of cosmic ray nuclei data (Maurin et al. 120011 12002J) . The same parameters as 
employed to study the antiproton flux (Barrau et aL I2002|l are used again in this analysis. 
The antideuterons density produced by evaporating pbhs per energy bin i\>f) obeys the 
following diffusion equation: 

dz * \dz 2 \d Z/ 

where q p " m (r, z, E) corresponds to the source term discussed at the end of this section. 
The total collision rate is given by T q = nHcr BH VQ where fJQ H is the total antideuteron 
cross-section with protons (Hagiwara et al. 120020 . The hydrogen density, assumed to be 
constant all over the disk, has been fixed to tih = 1 cm -3 . 

Performing Bessel transforms, all the quantities can be expanded over the orthogonal 
set of Bessel functions of zeroth order: 

oo 

i'D = T, N ^ P " m MC(x)) (9) 

i=l 

and the solution of the equation (jSJ) for antideuterons can be written as: 

AT D,prim tn \ ( ~~ V C L\ Ui\L) , . 

N * (0)=eXP hiTj ^smh( W 2) (10) 
where 

Vi {L) - 2J Q L cxp (&(L - *')) sinh (f (L z')) q Y m \z')dz> 

(V 2 c 2 1 1/2 

A l = 2hT l g e + V c + A-^cothj^ 



V c—„- K \ — \ r — ) ) \i>D(r,z,E) + 2h5(z)r B iP B (r > 0,E) = qi> ri ™(r,z,E) (8) 
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Energy changes (predominantly ionization losses, adiabatic losses and diffusive 
reacceleration) arc taken into account via a second order differential equation for 
N D, pnm ^ e g Eq ^ in Barrau et aL ^003, or Secs.3.6.1, 3.6.2 and 3.6.3 in Maurin 
et al. 1)2001(1 for further details). At variance with antiproton studies, in a first approx- 
imation, wc discarded the so-called tertiary term (corresponding to nonannihilating 
inelastic reaction, as given in Sec. 4 from Donato et al. (2000) which should be unimpor- 
tant at the considered energies since the binding energy of this nucleus is of about 2 MeV. 



The spatial distribution of pbh is a priori unknown. However, as these objects should 
have formed in the very early stages of the history of the Universe, it seems reasonable to 
assume that their distribution should be rather homogeneous. When the cosmic structures 
have formed, they should have followed the cold dark matter particles and we assume 
that they currently have the same distribution. As a consequence, the following profile 
for the PBHs distribution has been used (normalized to the local density): 

/('."> -i^^? (") 

where the core radius R c has been fixed to 3.5 kpc and Rq=8 kpc. This profile corre- 
sponds to the isothermal case with a spherical symmetry, the uncertainties on R c and 
the consequences of a possible flatness have been shown to be irrelevant in (Barrau et al. 

mm. 



4. Top of the atmosphere spectrum and experimental detection 

The flux is then solar modulated in the force field approximation with $ = 500 MV - 
corresponding to the solar minimum - and shown on Fig. [3 for a reasonable (po=160 
MeV/c, L=3 kpc) set of parameters at the top of atmosphere (TOA). The lower curve 
is the antideuteron background due to interactions of cosmic rays on the interstellar 
medium as given in (Donato et al. I2000[) whereas the upper curve is due to evaporating 
PBHs with a local density of 10 -33 g.cm~ 3 (allowed by the currently available upper 
limits (Barrau et al. I2002JI '!. Secondaries have been obtained in a two-zone diffusion 
model, with some simplifications: no convection and no energy losses have been included. 
However, as in the case of antiprotons, their effect should be marginal, while they are 
very important for primary fluxes, and the conclusions of the present analysis should not 
be substantially modified. To see all the computation steps, we refer the interested reader 
to Donato et al. H'2000| l: the procedure is basically the same as in this work, except for 
the production cross-sections that are simply deduced from the antiproton production 
cross-sections within a coalescence model with a fixed momentum (taken as 58 MeV, 
which corresponds to 116 MeV in our notation) instead of being computed by a Monte- 
Carlo method. The fundamental point is that this background becomes extremely small 
below a few GeV/n for kinematical reasons: the threshold for an antideuteron production 
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Fig. 2. TOA antidcutcron flux at solar minimum. The upper curve (left part) is from a 
pbh distribution with a local density of 10 -33 g.cm -3 and the lower curve (taken from 
Donato et al. I2000fl is from secondary processes. 

is E = 17 m p (total energy) in the laboratory 2.4 times higher than for antiproton 
production. The center of mass is, therefore, moving fast and it is very unlikely to 
produce an antideuteron at rest in the laboratory. It should be noted that the secondary 
antidcutcron background is only presented here to give a crude estimate of the expected 
"physical" background. In a forthcoming paper, we expect to study this secondary flux in 
much more detail, taking special care in the treatment of diffusion and the cross-sections. 

The number of events expected in the AMS experiment (Barrau 12001(1 onboard 
the International Space Station can be estimated, following Donato et al. I(20()0|l . 
Taking into account the geomagnetic rigidity cut-off below which the cosmic-ray flux 
is suppressed (as a function of the orbit parameters ), the acceptance of the detector 
and convoluting with the TOA spectrum, we obtain 7 events in three years between 500 
MeV/n and 2.7 GeV/n for the previously-given PBH density and the previously-given 
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Fig. 3. Parameter space (halo thickness L: 1-15 kpc; coalescence momentum pq: 60-285 
MeV/c; pbh density p & : 10~ 35 — 10~ 31 g.cm~ 3 ) within the AMS sensitivity (3 years of 
data). The allowed region lies below the surface. 

typical astrophysical and nuclear parameters. This is a quite low value which would be 
difficult to measure due to the possible mis-reconstruction of p or D events. Nevertheless, 
it should be emphasized that the situation is very different to that of antiprotons, as the 
limit here is not due to the unavoidable physical background but just to the instrument 
capability. Many uncertainties are still unremoved and can affect the primary flux more 
significantly than the secondary one. 

In order to be more quantitative, we performed a multi-variable analysis. Our model 
has a large set of free parameters: the astrophysical quantities used for propagation 
{Kq, 5, L,V c ,V a ), the local density p Q of pbhs and the nuclear coalescence momentum po 
for the formation of antideuterons. To evaluate the possible detection of a signal we chose 
the following strategy: as the main uncertainty from astrophysical processes comes from 
the halo thickness L, the other parameters were fixed to the value giving the smallest 
flux. This sub-set of parameters depends slightly on L and was varied as a function of 
L to ensure that whatever the thickness chosen the real minimum is reached. All the 
results are therefore conservative. The remaining variables p© , L and po are then varied 
within their allowed physical ranges: L between 1 and 15 kpc (see Barrau et al. H2002|l 
for the details), po between 60 and 280 MeV/c (depending on the experiments) and p Q 
on the largest scale matching the related experimental sensitivity. Two experiments were 
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Fig. 4. Parameter space (halo thickness L: 1-15 kpc; coalescence momentum po: 60-285 
MeV/c; PBH density p Q : 10~ 35 — 10~ 31 g.cm~ 3 ) within the GAPS sensitivity (3 years of 
data taking). The allowed region lies below the surface. 

investigated: the large spectrometer AMS (Barrau I2001JI which will take data over 3 
years from 2005 and the GAPS project (Mori et al. I2002f) . based on a clever design using 
X-ray desexcitation of exotic atoms. The allowed parameter space is given in FigJSJand 
FigEl the values of L, po and pq that can be explored by the considered experiment, 
without taking into account possible mis-reconstructions, are located below the surface. 
The sensitivity of AMS was taken to be 5.7 x 10 -8 m _2 sr _1 GeV/n _1 sec _1 between 500 
MeV/n and 2.7 GeV/n for three years of observations whereas the one of GAPS was 
chosen as 2.6 x 10~ 9 m _2 sr _1 GeV/n _1 sec _1 between 0.1 GeV/n and 0.4 GeV/n for the 
same duration (Mori et al. 12002)). To make the results easier to read, FigEJand FigElgfve 
the accessible densities of pbhs for AMS and GAPS with a fixed L (at the more reasonable 
value around 3 kpc) or a fixed po (at the more favoured value around 160 MeV/c). As 
expected, the primary flux is increasing linearly with the PBH density (at variance with 
the search for supcrsymmctric particles related to the square of p@, as a collision is 
involved) , linearly with the magnetic halo thickness (as the core radius R c is of the same 
order as L) and with the third power of the coalescence momentum (as the probability 
to create an antideuteron is related to a volume element in this space). The smallest 
detectable density of PBHs for the employed astrophysical and nuclear parameters is 
p & w 10~ 33 - 60 w 2.6 x 10~ 34 g.cm- 3 for AMS and p Q w lO" 34 ' 86 w 1.4 x 10~ 35 g.cm" 3 
for GAPS. It is much less than the best upper limit available nowadays p Q < 1.7 x 10 -33 



12 



Barrau et al.: Antideuterons from primordial black holes 



? E -32.25 
° -32.5 

Cn 

^-32.75 
'I -33 
x)-33. 25 
m -33.5 
^-33.75 
5 -34 
-34.25 
-34.5 




10 12 14 

halo size (kpc) 



. -32 
-32.25 
, -32.5 
-32.75 
-33 
-33.25 

-33.5 
-33.75 
-34 
-34.25 

-34.5 



50 



75 



1 00 



125 



I 50 



175 200 225 250 275 300 
coalescence momentum (MeV/c) 



Fig. 5. Upper plot: parameter space (pbh density vs halo thickness) within the AMS 
sensitivity for a fixed value of the coalescence momentum po = 160 MeV/c. Lower plot: 
parameter space (pbh density vs coalescence momentum) for a fixed value of the halo 
thickness L = 3 kpc. 

g.cm -3 and it should open an interesting window for discovery in the forthcoming years. 
If no antidcutcron is found, the upper limits will be significantly decreased, allowing 
stringent constraints on the spectrum of fluctuations in the Universe on very small scales. 
It should also be mentioned that, in spite of its much smaller acceptance, the PAMELA 
experiment (Adriani et al. I2002J) could supply interesting additional information thanks 
to its very low energy threshold, around 50 MeV/n. 

5. Discussion 

As recently pointed out in Donato et al. H2000J) . antideuterons seem to be a more promis- 
ing probe to look for exotic sources than antiprotons. In this preliminary study, we show 
that this should also be the case for PBHs, so that antideuterons may be the only probe 
to look for such objects. They should allow a great improvement in sensitivity during the 
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Fig. 6. Upper plot: parameter space (pbh density vs halo thickness) within the GAPS 
sensitivity for a fixed value of the coalescence momentum po = 160 MeV/c. Lower plot: 
parameter space (pbh density vs coalescence momentum) for a fixed value of the halo 
thickness L = 3 kpc. 

forthcoming years: a factor 6 better than the current upper limit for AMS and a factor 
of 40 for GAPS. 

Among the possible uncertainties mentioned in Barrau et al. (|2002|l . the most impor- 
tant one was, by far, the possible existence of a QCD halo around pbhs fHccklcr ll997fl . 
The latest studies seem to show that this effect should be much weaker (Mac Gibbon et 
at, in preparation) than expected in Cline et al. <|1999fl . The results given in this work 
should, therefore, be reliable from this point of view. 

Nevertheless, two points could make this picture a bit less exciting and deserve 
detailed studies. The first one is related to secondary antideuterons: the cross-sections 
used in this work could be slightly underestimated and some other processes could have 
to be taken into account (Protassov et al., in preparation). This could increase the 
background which should be considered with the same propagation model. The second 
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one is that the signal is extremely close to the one obtained with the annihilation of 
supcrsymmctric particles as the shape of the spectrum is mostly due to fragmentation 
processes. In the case of detection, it would be very difficult to distinguish between the 
two possible phenomena, unless collider data or indirect or direct neutralino dark matter 
searches have given enough information to fix the supcrsymmctric parameters. 

Acknowledgments. Wc would like to thank K. Protassov and R. Dupcrray for very 
interesting discussions about antideuteron cross-sections and C. Renault for her great 
help. 
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